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The impurity Green function and dynamical susceptibilties for the two-channel Anderson impurity 
model are calculated. An exact expression for the self-energy of the impurity Green function is 
derived. The imaginary part of the self-energy scales as ^J\uj/Tk\ for T ^ serving as a hallmark 
for non- Fermi behavior. The many-body resonance is pinned to a universal value l/(27rA) at = 0. 
Its shape becomes increasingly more symmetric for the Kondo-regimes of the model. The dynamical 
susceptibilities are governed by two energy scales Tk and Th and approach a constant value for 

— > 0, whereas relation x"('^) oc lo holds for the single channel model. 

PACS numbers: 71.10.Hf, 71.27.+a, 71.20.Eh 
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Introduction: In the last decade, the paradigm of a 
Fermi-liquid ground state in Heavy Fermion (HF) mate- 
rials 0, which originally gave these mainly Ce and Ura- 
nium base compounds their name, has become strongly 
questioned. In addition to magnetic and superconduct- 
ing phase transitions, many of these alloys show de- 
viations from the Fermi-liquid behavior in the normal 
phase. Typically, logarithmic behavior and power laws 
with anomalous exponents have been reported in the 
low-temperature thermodynamic and transport proper- 
ties . These unusual properties are attributed to charge 
and spin fluctuations on the partially filled /-shells. Cer- 
tain ballistic metallic point contacts, two-level systems 
and coupled quantum-dot have been argued to display 
similarly unusual properties. 

While the 4/-shell of Ce ions is usually occupied by 
one electron, the configuration of the 5/-shell of U ions 
fiuctuates between degenerate 5/^ and 5/^ crystal field 
states. Hund's rules coupling and crystal field splitting in 
a cubic symmetry yield a Kramer's degenerate Fg ground 
state |(t) for the 5/"^ and a quadrupolar S/^-Fa doublet 
la). On symmetry grounds, these states only hybridized 
with Fg conduction electron labelled by a spin tr and 
channel index The Hamiltonian of this two-channel 
Anderson impurity model (TC-SIAM) is given by 



kaa <y 



kacr 

where the operator cl,^^ creates an electron with energy 
Efcao-, momentum k, quadrupolar channel a = ± and 
spin a. Xm,m' = \m){m'\ is the usual Hubbard operator 
which destroys the ionic state |ra') and creates the ionic 
state |m). The energy Ecr is assigned to the magnetic 
doublet, the energy Ea to the quadrupolar doublet con- 
taining one electron less than \ a). The last term describes 
the hybridization between the /-shell and the conduction 
electron host. The external magnetic and channel fields. 



hg and he, couple to the local spin or channel flavor op- 
erator, Sz = J2a '^^•yal'i-, and = Yl,a o^^aa/'^ and are 
used to calculate the local spin- or channel susceptibil- 
ity = -TT = respectively, for hs,hc 0. 
This model was proposed in the context of UBeis by 
CoxHli for a novel and unusual Kondo effect: the non- 
magnetic S/^-Fa ground state undergoes a quadrupolar 
charge rather than a magnetic Kondo-screening. Using a 
Monte-Carlo approachjg, it was shown that (i) the mag- 
netic and quadrupolar moments are screened for T — > 0, 
(ii) this screening is described by two different energy 
scales Tfi > Ti ^ Tk, (iii) both channel and spin suscep- 
tibility behave as Xs/c oc lii(T) in the universality regime 
T <^Ti. Recently, this was confirmed by an exact Bethe 
ansatz solution[3 of the model It was also possi- 
ble to calculate the threshold exponents of the resolvent 
Im{z) = {M,Q\l/{z- H)\M,0) K Io{z - Eg)-""'" of the 
time evolution of the local levels M = a,a relative to the 
Fermi-see of the free electron gasQ; Eg is the threshold 
energy. In contrary to the non-Crossing Approximation 

(NCA)[illli[i3 predicting aM = 1/2, the true expo- 
nents are dependent on the occupation of Uc ~ '^^{Xcra) ■ 
Since the Green function Gaaiz) =<^ Xa^a\Xa.a ^ is 
obtained by a convolution of Icr{z) and Ia{z) in NCA 
[3, its spectral function would be given by paa{^) — 
- - a„) cx |c^r5"=(^-i) for 
— + using the exact threshold exponents [llj; B is 
the beta function. Hence, the spectral function should 
diverge for ric ^ 0, 1, indicating that vertex corrections 
are as important as the exact exponents om for obtaining 
the correct spectral function from /a/f fz)|l2|. 

In this Letter, the first accurate solution for the dy- 
namical properties of the model using Wilson's numerical 
renormalization-group method (NRG)^^ is presented. 
It is shown that (i) the scaling behavior of Paa{'^) is given 
by (1 — a± ^J\uj\/Tk) with a very weak ric dependence 
of the constants a±, (ii) the many-body resonance be- 
comes increasingly more symmetric in the Kondo regime 
in contrast to the NCA result, (iii) the total self-energy 
consists of a resonant level term and a local non-Fermi 



liquid self-energy with vanishing imaginary part at w = 0, 
(iv) the dynamical spin and channel susceptibility x"i^) 
approach a constant for < Tx and T <^ Tx, corre- 
sponding to logarithmically divergent x{T) for T — > 0. 

The Hamiltonian (jTJ is solved using the NRG[i^. The 
core of the NRG approach is a logarithmic energy dis- 
cretization of the conduction band around the Fermi en- 
ergy, uj^ = ±Z?A~" and A > 1, and a unitary transforma- 
tion of the base c ± onto a base such that the Hamil- 
tonian becomes tridiagonal. The first operator /oao- is 
defined as foaa = i E„(c^+„a + c<^-aa) M Only 
/oqo- couples directly to the impurity degrees of freedom. 
Eqn. is recasted as a double limit of a sequence of 
dimensionless NRG Hamiltonians: 



n= lim lim (DaA-^^-i^/^-^tvI, 

A^l+ W^oo L J 

with Da equal to D{1 + A"^)/2, and 



(2) 



Hn = A— 



Da 



N-l 



(7 (7Q. 

E E^'^^'- +H.C.} 

n— aa 

Va is related to the hybridization width DaVu = Va, 
AE = Eca ^ E„ - Ea, while the pre-factor K^^-^)/"^ 
guarantees that the low-lying excitations oi Hn arc of 
order one for all A^^3|- In the absence of a symmetry 
breaking channel or magnetic field, the Hamiltonian has 
three fixed points: the free orbital fixed point HpQ with 
AE = 0, y = 0, the two local moment fixed points H^j^j, 
\AE\ — ^ oo, \^ = 01 and an intermediate couphng 
or two-channel fixed point H^^ governing the universal- 
ity regime of the model for T ^ 0. HpQ is instable 
with respect to the two relevant operators X^^ and the 
hybridization, the local moment fixed point is instable 
with respect to the Kondo interaction Hk, which is ei- 
ther of spin type when AE < 0, or of channel type when 
AE > 0. There are two cross-over energy scales in the 
problem: the high temperature scale Th, associated with 
the screening of the upper doublet, and the low tem- 
perature scale Ti = Tk, the Kondo temperature, associ- 
ated with the screening of the lower doublet. The latter 
temperature is exponentially dependent on A/\Ea^\, and 
the Bethe ansatz predicts Tk = ^ exp[— ^^^^^ ] for 
the wide band limit. Introducing the hybridization func- 
tion Taaiz) = J2k ya\'^/{z-£kaa), the imaginary part at 
uj = —iS is denoted by A = Aaa — '^mVaai^i^)- 

Thermodynamic results: A symmetric and constant 
band density po = 1/{2D), a band width of £> = lOA, 
and a NRG discretization parameter A = 3.5 was chosen. 
All energies are measured in units of the Anderson width 
A. At each NRG iteration, Ng = 3000 states are kept. 
The impurity entropy Simp is defined as the difference of 
the entropy of ^ with or without the quantum impu- 
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FIG. 1: Entropy vs temperature for different values of AE = 
Ea<y. For |Ai5/A| <C 1, the Hamiltonian flows directly from 
the free orbital fixed point Hpp (Sfo ~ fcslog(4)) to the 
two channel fixed-point {Stc = log(2)/2). For large values 
|Ai?/A| ^ 1, the Hamiltonian crosses over to the local mo- 
ment fixed point {Slm ~ log(2)) first, before reaching the 
line of two channel fixed-points H^^. Tk associated with the 
quenching of the lower doublet is shown in the inset. Param- 
eters: D = 10, Ns = 3000, A = 3.5. 



values of \AE\ 



m the upper part of Fig. H and agrees 
excellently with the Bethe ansatz results [vj. For large 
values of |Ai?/A| 3> 1, the reduction of the entropy from 
log(4) — > log(2) occurs on the scale T/j, while the two- 
channel fixed point is approached below Tk 0, • Even 
though the two-channel fixed point is always character- 
ized by a residual entropy of i log(2), the NRG level spec- 
trum of the fixed-point is non-universal and dependent 
on AE or the occupation number ric- Since the model 
lacks particle-hole symmetry away from AE — 0, the ad- 
ditional marginal operator, H,n = K{J2aa flaofacro - 2) 
|l4| , leads to a line of fixed point Hamiltonians compris- 
ing of H*^{AE = 0) and Hjn- The constant K can be 
extracted from the NRG level spectra and contains the 
Uc-dependence Q|. 

The logarithmically divergent susceptibilities to the 
low temperature form x{T) ~ —^^^\og{T)+b was fitted 
to extract Tk- Tk can also be defined as the tempera- 
ture at which the effective moment of the lower doublet, 
plffiT) = rx(T), is reduced to /i^^/Tx) = 0.07 [3. 
These two temperatures coincide within the numerical 
error. The local spin susceptibility Xs for AE < 0, or 
channel susceptibility Xc for AE > is plotted in the 
lower graph of Fig. ^ for five different values of AE as 
a function of T/Tk, to illustrate the scaling. The inset 
shows the exponential behavior of Tk (NRG) and the ex- 
cellent agreement with the Bethe ansatz wide band limit 
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FIG. 2: NRG spectral functions for three different values of 
AE at T = in comparison with NCA spectral functions. 
The inset shows the behavior for < Tk- The dotted lines 
marks l/(27rA). For increasing |A_B|, the peak of the NRG 
spectrum becomes more symmetric and approaches a univer- 
sal curve. The NCA spectra are calculated at T = W'^Tk 
for AE = 0,3 and T = W'^Tk- The value of Tk are 
0.405285, 3.6 * 10"^3.3 * lO"'' for AE = 0, -3, -6. Parame- 
ters: D = 10, Ns = 3000, A = 2.2. 

cretization correction A\ — i log A for the Kondo 

coupling IQ: J Jeff = J/A\. 

Dynamical properties: Using the equation of mo- 
tion for Fermionic Green functions, z A\B ^= 
{{A,B})+ < [A, >, the commutators of the op- 
erators Xa^a- and Ck,aa with the Hamiltonian, [Xa^^,H] 
and [ck.aui H] respectively, and the local completness 
J2a "'^CTCT + J2a ^aa = 1 , it is Straight forward to derive 

\z — Eacr — Taa] ^Xa,a\Xa,a'^ = 

P„^+ <Oaa|^a,a> ■ (3) 

The expectation value of anti-commutator P^a = 
{Xa,a) + {Xa.a) is alway 1/2 in the absence of a sym- 
metry breaking field. The local composite operator Oaa- 

Oaa — ^ ^ a {X ^fj^uCk — a — a X —(j^ — Q-Ck — aa^ 

k 

)(4) 

k 

also transforms as Xa^^-- If the Green function Gqo-(z) 
is parameterized by the self-energy Tt^aiz), Gaa{z) — 
Paa[z - " aa ~ ^aa{z)]^^, tliis sclf-cnergy is given 
by the exact relation 

= Maa{z)/Gaa{z) , (5) 

defining Maa{z) =<^Oaa\Xcr,a^- In contrary to the 
single channel model'lSj, there exists no non-interacting 
limit, and the spectral weight is reduced to Paa = 1/2. 



It is straight forward to calculate the matrix elements 
of operator Oaa within the NRG. The spectral func- 
tions for Maa and Gaa are obtained from their Lehmann 
representation by broading the (5-function of the dis- 
cret spectrum on a logarithmic scale, i. e. 6{uj — E) 
g-6V4e-(iog('^/s)/6)=/(0f5|^|) ^ The broadening pa- 
rameter is usually chosen as 0.5 < 5 < 1, and we used 
b = 0.8. The self-energy T,{z) is then calculated through 
©■ Fig.[21shows the spectral function pacjoj) in compar- 
ison with the results obtain by the NCA|5|. The overall 
agreement between the NRG and NCA curves for the 
same parameter is reasonably good in spite of the incor- 
rect threshold exponents of the NCA. The NRG tends to 
overestimate the broadening of the high energy peaks at 
uj « AE, but produces highly accurate results for small 
frequencies \uj\ < Tk- The inset of Fig. |21 depicts the 
many-body resonance in the region for \u!/Tk\ < 10 on 
a rescaled axis \u;/Tk\^^'^- It highlights three major as- 
pects of spectral function: (i) p{lu) oc 1 — a±|w/Tif I^/^ for 
< Tk, (n) for increasing \ AE\, the spectral function 
becomes more symmetric around a; = and should ap- 
proach a universal curve, (iii) the peak value is pinned to 
l/(27rA) a.t u! — independent of AE. Even though the 
vertex correction to the NCA spectra vanishs only in the 
limit of large spin and channel degeneracy N and M^, 
the scaling behavior of the NCA spectra is apparently 
correct. The NCA peak, however, saturates at 7r/(16A) 
for r ^ in and remains asymmetric for |Ai5| ^ A. 

Fig. 13 shows the corresponding imaginary part of the 
NRG self-energy I]{u! — iS): it vanishes at u = 0, inde- 
pendent of AE, which is a neccessary condition for the 
pinning of the spectral function for all values of AE. In 
addition, 3fieE(0) + A + Eaa = 0, leading to the pinning 
of the resonance to = 0, depicted in Fig. |21 In the 
scaling regime \ll;/Tk\ < T/i, 3mE(w — iS) oc y^\uj/TK\. 
This region extends to multiples of \uj/Tk\ for Th ^ Tk. 
The definition of the self-energy via the parameteriza- 
tion of Gaaiz) turns out to be very useful for the un- 
derstanding of its properties. The conformal field the ory 
result for the T-matrix of two-channel Kondo model llTl 
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FIG. 4: Imaginary part of x"('^) {^) scaled to position 
and value of the finite frequency maximum Umax ~ 7)i 
as displayed in the inset; the values of x'si'-^max) are 
0.0126, 0.0046, 0.00235 for AE = 2, 4, 6. The scaling curve 
(b) for vs uj/Tk at T = for AE > 0. Parameters as 

Fig.H 



conveys to the TC-SIAM. The self-energy 'S{z) exhibits 
similar features as the correlation self-energy 'S^ (z) — 
E(z) — r(z) of the single-channel model: its imaginary 
part also vanishes for w ^ and its scaling behavior re- 
flects the nature of the local electronic excitation. While 
3mE(a; — i5) cx (w/Tr-)^ in the single channel model, 
the scaling as 3mE(w — iS) cx ■^/jZT/T/ff represents the 
hallmark of a local non-Fermi liquid. 

The imaginary parts of the dynamical impurity spin 
and quadrupolar or charge susceptibility, Xsi^) =^ 
Sz\sz^ and Xc(w) =<S^Tz\tz^ respectively, are plotted in 
Fig.Elfor T = and positive AE > 0. Since the Hamilto- 
nian is symmetric with respect to a particle-hole transfor- 
mation, and simmultanious AE ~^ — AE and interchang- 
ing the spin and quadrupolar indices, Fig.0]also contains 
the information for AE < 0. Both x" show two frequency 
regimes. In contrast to the single channel model where 
the relation x" oc uj holds, x's/ci^) approachs a constant 
at low frequencies, \uj\ < Tk- This corresponds to the 
logarithmic divergence of both x(T) for T < Tk, and the 
slope of the logarithm is given by this constant. A scal- 
ing curve is found for the charge susceptibility (or spin 
susceptibility for AE < 0) when plotted in dimensionless 
units of the low temperature scale Tk, as shown in Fig. 01 
(b). The scaling holds for oj/T^ < 1 above which non- 
universal band edge effects govern the physics. The inter- 
mediate frequency range Tk < < Th is characterised 
by a power law behavior: Xci'^) °^ (w/Tj^)"^-^. The 
peak position in X^w) scales with Th and corresponds to 
AE for large AE 6, 7]. The values of x'si^rnax) become 
increasing smaller with reduction of the occupance of the 



spin-doublet. 

Discussion and Conclusion: The dynamical proper- 
ties of the two-channel Anderson model were accuratly 
calculated for T — using the NRG approach. The 
analytical form of the spectral function is given by 
l/(27rA)(l - a^^ujjTK) for |w| < Tk- Its peak is 
pinned to a universal value, and the self-energy scales 
as ^Jijj/Tk independent of n^. Since the self-energy 
E(z) does not contain the resonant level part A^cn and 
3mE(a; — «i5) — *■ for w — *■ 0, its anomalous scaling serves 
as a hallmark for local non-Fermi liquid behaviour. The 
NCA result for the spectral function agrees remarkably 
well with the accurate NRG result. The two energy scales 
characterising the screening of the moments of upper and 
lower doublets also appear as charge exitation energy and 
as effective width of the many-body resonance at w = in 
the spectral function. At finite temperatures - not shown 
here - the many-body resonance decreases for inceasing 
temperature and vanishes for T ^ Tk- The dynamical 
susceptibility of the lower doublet maps to a universal 
scaling curve for lo < min{I?, TfJ\ and T = 0, while a 
peak at oJmax ~ Th is found in the dynamical suscepti- 
bility of the upper doublet. 
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